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In this paper we present some results on finite energy solutions in an SU(2) Yang-Mills (YM) 
and YM-Higgs (YMH) theory (with the Higgs field in the adjoint representation) coupled to a 
massive dilaton field. The present work is an extension of previous investigations with a massless 
dilaton field 0|. In Refs. [Q, it has been found that the SU(2) dilaton-YM (DYM) theory 
admits static, finite energy (‘particle-like’) solutions (absent in the pure YM case). They are in 
close analogy to the particle-like solutions found by Bartnik and McKinnon in an Einstein-YM 
(EYM) system [^. In Ref. |Q] it has been shown that in the SU(2) dilaton-YMH (DYMH) theory 
in addition to the analogue of the (nonabelian) ‘t Hooft-Polyakov monopole there is a discrete 
family of finite energy solutions, which can be interpreted as radial excitations of the monopole. 
The mass scale of the radial excitations is inversely proportional to the dilaton coupling. 

In Ref. Q it has also been found that there is a maximal dilaton coupling, Omax, above which 
only an abelian solution exists. Although the abelian solution (which has a simple analytical 
form) is singular at the origin, its total energy is finite. The numerical results of Ref. indicate 
that the nonabelian monopole merges with the abelian one for a critical value of the dilaton 
coupling, a. The critical dilaton coupling, Ocrit, depends on the value of the Higgs self-coupling 
strength, (5. As found in Ref. Q when (5 is sufficiently small the largest possible value of the 
dilaton coupling, Omax! for which a nonabelian solution still exists is different from the critical 
value, acrit, he. we have Omax > «crit- This implies that if a G [acritj ctmax] a bifurcation takes 
place and there are two different monopole solutions for the same value of a. All these findings 
are again in close analogy with the results of Refs. ii found in the EYMH case. 

A massless dilaton, which necessarily appears in string theories i I®, violates the equiv¬ 
alence principle, and therefore its (dimensionful) coupling strength is expected to be extremely 
weak, of the order of 1/Mpi where Mpi is the Planck mass. It is very natural to assume, however, 
that the dilaton gets a mass (possibly related to supersymmetry breaking) then, however, there 
is no strong experimental constraint on the dilaton coupling. Therefore it might be of some 
importance to study the effect of mass of the dilaton in DYMH theories. Regular and black hole 


solutions in EYM theories coupled to a massive dilaton and axion have investigated in Ref. [ IT |. 

We have carried out a rather thorough numerical investigation and we have found that even 
if the dilaton is massive, most phenoma associated with the presence of the dilaton found in 
the massless case persist (the existence of radial excitations, Omax / «crit)- Also we have good 
numerical evidence that Ocrit (where the solution become abelian) is independent of the dilaton 
mass. We have also investigated in detail the limit m —>■ oo of the mass of the dilaton and we 
were able to show its expected decoupling. Our numerical results show that Omax grows as m 
increases, consistently with the expectation that for m —> oo a ma x(m) —> oo. 

We have also found that not just a single maximal dilaton coupling, Omax exists, but there 
are several local extrema of a too, e.g. amaxi > Q ma xo > Ocrit > CKmim • These bifurcation points 
are more clearly distinguishable as the dilaton mass becomes larger. This implies that there are 
certain values of the dilaton coupling, a, where three or even more different monopole solutions 
exists for the same value of a with different masses. 

By using the minimal spherically symmetric and static ansatz 


Ag — 0 , A“ — Caik 


1 - W{r) 


= H{r) 


ip = if{r) , 


( 1 ) 


where A“ is the gauge potential, <I>“ is the Higgs triplet and ip is the dilaton field, the reduced 
energy functional reads: 


E(a, (3\m^) = Jdr (w'^ + 
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Varying the energy functional @ we obtain the field equations: 
{r\')' = 2ae2“^ (w'^ + 

(e^^^lV')' = W 


- 1 


+ H^ 


{r^H'y = 2H _ ^ 2 ^^ _ 


By introducing the following variables: 


R = re := rep, v := aip, f" = ^ = — 

dr p 


the field equations (^) can be written as: 

R = 1 — Ru, 




m 


- ^R\H^ - vy + ■yR^iye^'' , 




W + Wu = W 


- 1 


R2 


+ H^ 


{R^HY+R^Hu = hUw^ + -v^) 


The field equations (^) written in this form are distinguished by the fact that for m = 
contain only variables which are invariant under the dilatational symmetry 

rvr* 1 

r —> e r, (p —> P + c, 


(3a) 

(3b) 

(3c) 

( 4 ) 

(5a) 

(5b) 

(5c) 

(5d) 

0 they 

( 6 ) 


which is only broken by the presence of the mass term. 

Solutions regular at r = 0 (or equivalently at r = 0) must satisfy the following boundary 
conditions: 

R = ar + O(r^), IT = 1 — br^ + O(r^), 

R = t- yr^ + 0(r®) , i^ = i^o + + O^) , 

with ui = i 46^ — 


m 

- H-t'oe' 

12 y 3 


i.e. solutions with a regular origin depend on three free parameters (a, b, uq). The asymptotic 
behaviour (r —> 00 ) of regular solutions is almost identical to those of the m = 0 case (Eqs. (10) 
in Ref. |^), i.e. for r —> 00 the b.c. corresponding to regular solutions are 


O' 1 

W(r) = Be-^ (^l + 0(^)^ , 

R(r) = 1 ——e~^^ f 1 + 0(-)^ for /3 < 2 
r \ r J 

^(0 = 1 - (^ + ^( 7 )) for/3>2, 


(8a) 

(8b) 

(8c) 

(8d) 


where {B, C) are free parameters and we have set n = 1. We note that p{r —> 00 ) does not 
tend to zero as which is the expected asymptotic behaviour of a massive field. This 
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‘pathological’ behaviour and the missing free parameter are due to the fact that in Eq. (|3|a) as 
IE —> 0 the 0(l/r^) term dominates. (The asymptotic behaviour of the Higgs field for /3 > 2 in 
Eq. d^) is determined by the nonlinear terms in a similar way.) 

With the use of the field eqs. i) and the boundary conditions (|^ one can easily show that 
the energy of the solution, E(q:,/ 9^,m^), is a monotonically increasing function of and 
and a monotonically decreasing function of a. 

In the DYM limit of the theory, i.e. (u ^ 0 and H > 0) there are no monopoles, but our 
numerical results clearly indicate the existence of globally regular solutions with zero magnetic 
charge. The existence of these solutions is not surprising, as they are simply the generalizations 
of those found (also numerically) in the massless case (m = 0) in Refs. |]|, Our results for 
the m = 0 case are in complete agreement with those of Refs. |j^, Q]. 

We have found only a slight dependence on the dilaton mass. Although our numerical 
solutions exhibit clearly that the dilaton tends to zero more and more as m increases, W varies 
very little for the considered range of values of m ( 0 < m < 12) as it can be seen on Fig. [l|. 
In the m ^ oo limit no nontrivial regular solution survives, because the dilaton tends to zero 

W, 1/0 W, 1/0 




Figure 1: dependence of the first two excited solution in the DYM limit. 


everywhere forcing W = 1. This limit is, however, singular and we expect that for finite values 
of m (no matter how large) an infinity family of regular solutions still exists. 

Before considering nonabelian monopole solutions, we discuss first the abelian case, i.e. 
= (^ab) W = 0, H = V {a Dirac monopole coupled to a massive dilaton field). There is a finite 
energy solution, with a (logarithmic) singularity at the origin. In the massless case (/^ab = 
and it is given by 




(0) = -iln 


l + “ 


( 9 ) 


a 


r 


while for m 7 ^ 0 its analytic form is not known. Note that for abelian solutions the only relevant 
parameter is ma (for m 7 ^ 0, a 7 ^ 0). It seems very likely that the abelian solution is unique for 
any value of m. 

Near the origin the power series expansion of ip^h can be written as 


, m‘^ 9 s, r T 9 5 9s 9 s s 

aipah = (1 + —rr) In - + ar + (-a^ - —m^)r^ + ... ( 10 ) 

4 a 2 lb 

where o is a free parameter. So clearly if r <C 1/m the mass term in ( 0 ) is negligible and the 
dilaton is well approximated by the regular solution with m = 0, see Fig. 

Next we investigate the nonabelian solutions. In the asymptotic region, r —> 00 , the be¬ 
haviour of the nonabelian solutions is determined by linearizing the field Eqs. (P) around the 
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Figure 2: w? dependence of the abelian solution for a = 1. See Table for their parameters. 


Table 1: Parameters of the abelian solution for a = 1 



a 

E 

0 

1.000000 

1.000000 

1 

0.055585 

1.304817 

10 

-2.164717 

1.724591 

100 

-10.124077 

2.594755 

1000 

-38.701493 

4.278556 


critical point {W = 0, = u, (^ = 0} for r —> oo, yielding Eqs. (^). As in this case ^9 —> 0 

Rk. r, r K, T, W and [H — l)r are exponential functions of r. 

We solved the field equations @ numerically by a combination of a Runge-Kutta (RK) 
method (4th order, adaptive stepsize) integrating out from r = 0 to some tq and solving by 
iteration the pertinent system of integral equations from r = oo to tq. In order to desingularize 
the field equations (D at the origin resp. at r = oo, we have introduced the following set of 
variables: {{W -l)/R, W, H, (rH)', 0, ^,R,} resp. {W± = W±W,h± = f3h ± H+, 
= mRif ± y, R — t} where h = {H — l)R, H+ = (rH)' — 1, Y = ip + Rep. 

(For /3 = 0 we used R^H instead of /i_, and for m = 0 we replaced p- by R?ip.) The corre¬ 
sponding eqs. for (3^0 and m ^ 0 are the following: 


VF± =T 


/i± ==F 


A/t 
l/r± 

H/t 
1/t± 


w[-^Y-{}ilRY‘l)yaWp 






R 


2R 






"l/r 


= =F 


/i/t± 


a 


2W^ + 


(W^ - 1)2 /32/i2 


(h/R + 2)^ 


R-t = 


4 

+m^K{e‘^^^ - 1) T mKap] 

l/r 
l/r- 


rl/r 

/ aRipR dt + Cr-j 

J I/t— 


(lla) 


(lib) 


( 11 c) 

(lid) 


where r+ = oo, r_ = tq. The C^s are arbitrary constants determined by the boundary conditions 
( 0 , 1 ) ensuring regularity at r = 0 and at r = oo. In order to supress the divergent modes in 


5 

























{VF+,/i+, (^+1 at r = oo one has to choose Cw+:Ch^-,C^p^ = 0. The remaining four constants 
{Cw-,Ch_,C^_,CR-r} are determined by Cw. = W-{to), Ch_ = h-{To), C^_ = v^-(to), 
Cr-t = R{to) — To, where {W-{to), h-{To), (/?_(ro), R{to) — tq} are given by the RK procedure, 
guaranteeing regularity of the functions at r = 0. The integral eqs. (^) are then solved by 
iteration yielding a regular solution for r G [tq , oo]. Then the shooting parameters for the RK 
are adjusted so that the resulting values {lT+(ro), /i+(to), (^+(ro)} match those calculated from 
the solution of the integral equations. 

Let us now present our numerical results in some detail. First of all the fundamental non- 
abelian monopole solutions, which tend to the’t Hooft-Polyakov monopole as a ^ 0 exists only 
up to a maximal value of the dilaton coupling strength a, rn). For a > Omax only the 

abelian solution {ip = VF = 0, FT = 1} seems to exist. The a dependence of the solution 
is, however, quite complex as the nonabelian solution does not necessarily cease to exist for 
a = Omax- In fact that value of a (which value we call critical, Ocrit) for which the nonabelian 
solution merges with the abelian one, does not necessarily coincide with Omax- As a approaches 
its critical value, Ocritj the dilaton field becomes divergent at the origin, whereas IT(r) —> 0 and 
H{r) —> 1. This behaviour is shown on Fig. where the fundamental monopole is plotted as a 
varies from 0 to Ocrit- Note that although the solution becomes singular its total energy remains 
finite. 


W, H 1/0 






Figure 3: The a dependence of the fundamental monopole solution for m = 1, = 0.1. Above 

the exterior, below the interior solution is depicted. 


The critical value of a depends strongly on /3, but it does not seem to depend on the mass 
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Table 2: Parameters of the fundamental monopole solution for m = 1, = 0.1 


a 

a 

b 

-Po 

aE 

0 

0.4821127 

0.2381420 

0.0000000 

0.0000000 

1 

0.4481335 

0.2370648 

0.2787202 

1.0605540 

1.35 

0.4021605 

0.2368443 

0.6065505 

1.3627414 

1.4 

0.3874162 

0.2378367 

0.7807487 

1.3967137 

1.419 

0.3750471 

0.2400726 

1.0195001 

1.4078919 

1.419 

0.3719369 

0.2411243 

1.1136470 

1.4078993 

1.4 

0.3650065 

0.2464302 

1.5649554 

1.3999542 

1.35 

0.3662630 

0.2550531 

3.8698184 

1.3849301 

1.3495957 

0.3662688 

0.2546880 

9.3420476 

1.3848394 


of the dilaton, m. This latter phenomenon can be understood as follows. By performing the 
transformation (^) we define (^sh(?’) = — ^po where ipo := (^(0) guaranteeing (/?sh(0) = 0 and 

then the field eq. takes the form: 

(fVsh)' = 2ae2“^=‘' + + , (12) 

where f = and the prime denotes now derivation with respect to f. Note that v^sh(f^) 

—ipo as f ^ oo. The remaining field eqs. (^ 3 ,^ 0 ) are form invariant. 

Using the shifted field, (/?sh) in order to obtain a globally regular solution with some fixed 
asymptotic value of y?sh) -¥^ 0 ) one also has to vary a in addition to the shooting parameters a, 
b. (One still has to supress three divergent modes). Then the critical value of a is defined by 
taking the limit ipo —> — 00 . Then in Eq. (^^ for any, fixed value of f the terms proportional 
to m? tend to zero as ipo — 00 . From this it follows that Ocrit is determined by the massless 
equations. Note that it is essential to keep f fixed when taking the ipo —> —00 limit. In this limit 
one obtains nontrivial {W{f), H{f)} which are globally regular and a logarithmically divergent 
ip{f) as f —> 00 . Since for any fixed value of f r —> 0 as (^0 “OO this solution will be referred to 
as the ‘interior’ solution. It corresponds to a singular solution of Eqs. (P), IU(r) = 0, H{r) = 1 
for r 7 ^ 0 with a logarithmically diverging dilaton field at r = 0 which we call ‘exterior’ solution. 
Note that while the interior solution is nonabelian, the exterior one is, however, abelian. 


Table 3; and dependence of ccmax) «miin «c 



m2 

13^=0 

0.07 

0.1 

1 

O^max 

0 

0.1 

0.5 

1 

5 

1.408800 

1.411778 

1.471277 

1.586124 

2.314522 

1.361925 

1.361912 

1.374598 

1.439274 

1.965262 

1.359389 

1.419335 

1.924131 

1.206183 

1.206184 

1.209233 

1.244148 

1.612840 

C^min 

0 

0.1 

0.5 

1 

5 

1.399088 

1.39581 

1.361900 

1.361898 

1.361843 

1.361764 

1.36137 

1.349560 

1.349553 

1.34947 

1.20618034 

1.20617980 

0^max2 

0 

0.1 

0.5 

1 

5 

— 

1.361902156 

1.361902253 

1.349595705 

1.349596618 

— 


0 

0.1 

0.5 

1 

5 

1.3993823 

1.3993823 

1.3993823 

1.3993823 

1.3993823 

1.361902 

1.361902 

1.361902 

1.361902154 

1.361902154 

1.349595679 

1.349595679 

1.349595679 

1.349595679 

1.349595679 

1.20618036 

1.20618036 

1.20618036 

1.20618036 

1.20618036 


With our numerical procedure we found that in fact there is a local minimum of a, aminj 
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and that even a second local maximum of a exists. We remark that some values (for m = 0) 
differ from our previous results as given in Table 2 of ref. Q. The present data is more precise 
due to the use of the integral equations around infinity. 

The numerical data show that the local minimum of a exists for a wide range of jS, especially 
in the case of large dilaton masses. This local minimum exists even in the m = 0 case, but for 
larger m the difference between the different extrema of a and of Ocrit is more visible. This is 
the reason why amin has not been found in Ref. 0, where the present matching procedure has 
not been used. It is quite possible that even more local extrema of a exist but our numerical 
data is not precise enough to see them. To investigate this question a better method is needed 
as one has to approach Ocrit very closely where our numerical errors become too large. As (5 
increases (with m fixed) Omin and a;max 2 gets closer. The mass and the j3 dependence of the 
different extrema of a together with those of Ocrit are given in Table |3|. We also notice that amax 
grows quite rapidly as m increases while the critical alpha remains unchanged. This behaviour 



a 

1/00 



b 



E 



Figure 4: The a dependence of the shooting parameters a, 6, l/4>o and the energy of the 
fundamental monopole with (P = 0.1 for some m?. 

is illustrated on Fig. ^ where the a dependence of the shooting parameters a, b, and the 
energy of the solution, E, for different values of m are plotted. It is clearly visible that they 
all become more and more constant for a ^ Omax as m grows. So the solution becomes more 
and more independent of a for a G [0, Omax] as m increases and gets closer and closer to the’t 
Hooft-Polyakov monopole. The numerical results indicate that in the limit m —> oo the functions 
W and H approach those corresponding to the’t Hooft-Polyakov monopole while <y9 —> 0. 


8 
















To understand the m ^ oo limit somewhat better let us define <f{r) by (^ = a(p/m?. As¬ 
suming now that ip{r) ^ 1 and keeping only the leading terms in vn? it follows from the field 
equations that W « Wthp, H k, i^thp where Wthp, -f^thp denote the ’t Hooft-Polyakov solution 
and (p K (^inf where 

<^inf = + (Wl^ - l)Vr2 - - 1)^^ . 

The numerical results show that already for m? = 10 (p^ai approximates (p very well, see Fig. |^. 
The shooting parameters also tend to the appropriate (’t Hooft-Polyakov) values. Fig. ^ also 
exhibits the tendency that the shooting parameters and the energy of the solution depends less 
and less on a as m increases. 

So the above result supports the expectation that in the limit m ^ oo the dilaton decouples 
while tending to zero and W, H converging to the’t Hooft-Polyakov monopole solution. Note 
that the growing of the maximal a, Omaxj as m becomes larger is also consistent with the 
expected decoupling. 

0 

- 0.1 
- 0.2 
- 0.3 
- 0.4 
- 0.5 
- 0.6 
- 0.7 

0 0.2 0.4 0.6 0.8 1 

r! (r -|- 1) 

Figure 5: The m? depedence of the fundamental monopole solution for a = 1, fp' = 0.1. 
(p = rr?p>la and ’tHP denotes the ’t Hooft-Polyakov monopole solution. The corresponding 
parameters are given in Table ^ 



W, H 



Table 4: Parameters of the fundamental monopole solution for a = 1, = 0.1, together with 

the corresponding values for the’t Hooft-Polyakov monopole. The starred value stands for the 
limiting value of —(^(0) for m? ^ oo i.e. for —<^inf (0) = 126^ — /3^/4. 



a 

b 

-po 

E 

0 

0.4292657 

0.2386866 

0.8816368 

0.9001767 

0.1 

0.4321814 

0.2384289 

0.6263968 

0.9845152 

1 

0.4481335 

0.2370648 

0.2787202 

1.0605536 

10 

0.4726173 

0.2373279 

0.0539949 

1.0988876 

100 

0.4808578 

0.2380086 

0.0063885 

1.1054006 

’tHP 

0.4821127 

0.2381420 

0.6555393* 

1.1061955 


In addition to the fundamental monopole there is a discrete family of nonabelian solutions 
characterized by the number of zeros of W, n. They can be interpreted as radial excitations of 
the fundamental (n = 0) monopole. See Fig. One can understand these oscillating solutions 
by linearizing the field eqs. @ around {W = 0, H = v, (p = l/aln(r/a)} i.e. {R ^ a) in the 
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W, H 



ln(r) 


0 
-5 
-10 
-15 
-20 
-25 
-30 
-35 
-40 
-45 

-30 -25 -20 -15 -10 -5 0 5 10 

ln(r) 



Figure 6: The first five excited monopole solution for a = 0.5, = 0.1, w? = 1. Their 

parameters are given in Table 


Table 5: The parameters of the fundamental and first five excited monopole solutions for a = 0.5, 
/32 = 0.1, m = 1 


n 

a 

b 

-Po 

aE 

0 

0.4746849 

0.2378870 

0.1120060 

0.5478783 

1 

0.8255632 

1.5068349 

7.2301508 

1.1638519 

2 

0.8669493 

1.5525918 

16.3566790 

1.1741097 

3 

0.8676030 

1.5530864 

25.2510967 

1.1742216 

4 

0.8676121 

1.5530923 

34.1370513 

1.1742229 

5 

0.8676122 

1.5530924 

43.0228028 

1.1742229 


region r <^llm and r a corresponding to r ~ alnr+const. (For m = 0 {W = 0, H = v, 
/? = a, z/ = 1/a} is a critical point of Eqs. (|5|).) The linearization yields the following behaviour: 
{W{t) oc e'^'^'^, {H{t) — 1)t oc oc e'^'^’’}; where the exponents. A*, are given as 

X]y = l/a(—1/2 ± \/ot^ ~ 3/4), 

! - ( 13 ) 

\h = l/a(l/2 ± Y'a2/32 + 1/4), = l/a(l, -2). 

The amplitude of the oscillations of W tends to zero as 1/y/r and their frequency is independent 
from the mass of the dilaton. From Eqs. ( p!^ one expects that solutions with an oscillating 
W and a logarithmically growing dilaton (in the region where r <Cmin{a, 1/m}) exist only for 
< 3/4 . For larger values of r the behaviour of ip changes from p oc ln(r) to y? oc 1/r^ and 
W stops to oscillate. We expect that globally regular solutions with an arbitrary number of 
oscillations exist for a^ < 3/4. 

In the massless (m = 0) case we have found excited solutions (with n 7 ^ 0) only up to 
a = independently of (5 and n. But with an increasing mass we have also observed 

solutions with n = 1 for a > \/3/2 (E.g. a = 0.98, m^ = 10 for = 0.1). In fact there is also a 
maximal and a critical value of a for the excitations similarly to the n = 0 case. Our numerical 
data is, however, not precise enough to allow for a good determination of acritj (for n 7 ^ 0 ) but 
they are consistent with acrit = \/3/2 independently of m and n. This is in agreement with the 
argument given for the mass independence of acrit- We remark that the maximal value of a for 
solutions with n = 1 has always been less than one. The energy of the excitations increases 
with n and it is always larger than that of the fundamental monopole. The first zero of an 
oscillating solution moves closer and closer to r = 0 as the number of oscillations increases. The 


10 






























excited solutions depend only slightly on the dilaton mass, similarly to what has been seen in 
the massive DYM theory (u —> 0 limit). 

We expect that the excited solutions are unstable against spherical perturbations, because 
of the established instability of their counterparts in the massless DYM theory [Q, We also 
expect the fundamental monopole solutions to become unstable after the bifurcation point at 
c, analogously to the case of gravitating monopoles |^|. 


a = On 
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